The Edge of the Milky Way Stellar Disk Revealed Using Clump 
Giant Stars as Distance Indicators 



D. R. K. Saito^ J. Alonso-Garcia^ P. W. Lucas"^ and M. Hempel^ 

^ Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, 
Vicuna Mackenna 4860, Casilla 306, Santiago 22, Chile 
^ Vatican Observatory, Vatican City State V-00120, Italy 
^European Southern Observatory, Vitacura 3107, Santiago, Chile 
^ Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield 

ALIO 9AB, UK 



ABSTRACT 



We use the clump giants of the disk as standard candles calibrated from 
Hipparcos parallaxes in order to map their distribution with two new near-IR 
surveys of the Galactic plane: UKIDSS-GPS and VVV. We explore different 
selection cuts of clump giants. We conclude that there is an edge of the stellar 
disk of the Milky Way at i? = 13.9 ±0.5 kpc along various lines of sight across the 
galaxy. The effect of the warp is considered, taking fields at different longitudes 
and above and below the plane. We demonstrate that the edge of the stellar 
disk of the Milky Way can now be mapped in the near infrared in order to test 
different models, and to establish our own place within the galaxy. 

Subject headings: Galaxy: disk — Galaxy: structure — Stars: distances — Stars: 
late- type 



1. Introduction 



Most spiral galaxies have rather sharp edges in their stellar disks (Ivan der Kruitlll979 



20071 ). The situation for the Milk y Way is unclea r. The disk of our galaxy is well represented 
by an exponential density profile ( iFreemarul 19701 ) . with s ome evidence for a stellar disk cut-off 



at a radius of ~ 14 kpc towards the anticentre regio n (IRobin et al. 



for a warped and flaring stellar disk out to 23 kpc (IMomany et al 



19921), but other claims 



20061 ). Based upon the 



scaling relation established by other galaxies, the radial cut off should be between 10 and 
25 kpc ( iPohlen et al.l 120041 : iKregel fc van der KruitI 120041 ). Here we report observations of 
"clump red giant" stars, which are good distance indicators (lAlvesll2000l : iPaczynski &: Stanek 
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19981 ). along 10 lines of sight in the disk of the Milky Way. The data reveal an edge to the 
disk at -R = 13.9 ± 0.6 kpc. Why there should be a sharp edge to the stars, while the ga s 
profile is much more extended and does not show such an edge (IKalberla &: KerpI |2009| ) , 
rema ins unclear, but a critical test would be the measurement of the gas and star density in 
situ kennicuttlEgsgl l. 



There are two new deep surveys of the Milky Way disk in the near-infrared: the UKIDSS 
Galactic Plane Survey (GPS), that is mapping the Northern disk ( iLucas et al.ll2008l ). and the 
VISTA Variables in the Via Ldctea (VV V) Survey, that is rn apping 520 square degrees in the 
Southern disk and bulge of our Galaxy (IMinniti et al.ll2010[ ). While previous observations in 
the Galactic plane have been limited by crowding, source dimness, and interstellar extinction, 
these new public surveys allow us to pierce through the whole disk, reaching red giants in 
the horizontal branch (clump giants) located at the other side of our galaxy for the first 
time in different near-IR pass-bands. The VVV survey observes in the Z (0.87 microns), Y 
(1.02 microns), J (1.25 microns), H (1.64 microns), and Kg (2.14 microns) bands, while the 
UKIDSS-GPS is observing in the J (1.25 microns), H (1.64 microns), and K (2.14 microns) 



band s; with the photometry uniformly calibrated in the 2MASS system (iSkrutskie et al. 
2006h . 



Clump giants are excellent tools to study the structure of our galaxy because of two 
main reasons, discussed below: (i) they are good distance indicators, and (ii) they are good 
tracers of the mass in the disks of galaxies. The UKIDSS-GPS and VVV Surveys give the 
largest homogeneous census of clump gi ants out to well beyo nd the extent of the Milky 
Way stellar disk. The limiting DoPhot (jSchechter et al.l Il993l ) magnitude of single epoch 
VVV images processed by the Cambridge Astronomical Survey Unit (CASU) pipeline vl.O 
is Ks = 17.5 in the disk fields, and the DoPhot photometry on the combined epochs allow 
to reach approximately J = 20 mag, and Ks = 18.5 mag (Fig. [T]), slightly fainter than the 
UKIDSS-GPS, and much fai nter than the near-IR surveys 2MASS and DENIS, for which 
the limit is Ks = 14.3 mag (ISkrutskie et al.ll2006l ). The distance probed along the line of 
sight depends on the reddening of the fields. For example, in disk fields with low extinction 
{Ay < 3 mag) the UKIDSS-GPS and VVV surveys would see clump giants beyond 50 kpc. 
Therefore we can search for the edge of the disk of our galaxy using well calibrated standard 
candles. 



2. Observational Data 



The clump giants with known parallaxe s measured by the Hippar cos satellite are the 
best collectively calibrated standard candles (jPaczynski fc Staneklll998l ). 
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Fig. 1. — Selection of giant stars used in the distance calculation. The left-hand panel shows 
the colour-magnitude diagram (CMD) of VVV field d003 as a density map. The adopted 
color-cut is shown, which is similar to the cut for other fields, even though the location of 
the edge feature is relatively insensitive to these cuts. The CMD in the central panel and 
the colour-colour diagram in the top-right panel show in dark gray the selection of giant 
stars used in the distance calculation. Dashed lines in the CMD mark the distance in kpc 
for the selected stars according with equation 2. The bottom-right panel shows the distance 
distribution and respective derivative curve, with the outermost the minimum derivative 
marked. 



Figure [T] shows a colour-colour and colour- diagram magnitude for the VVV field d003. 
In order to select the giant stars a colour cut was applied in the Kg vs. {J — Kg), as shown in 
left-hand and central panels of Fig. [H Not all selected stars are clump giants, but the pres- 
ence of sub-giant branch and red-giant branch stars, or even main-sequence stars at fainter 
magnitudes, contribute to the smooth underlying background. These different population do 
not affect the location of any sharp edges of the stellar distribution, which can however be 
well defined by the clump giants. The effect of reddening in the selection of clump giants was 
also tested using different selection criteria both in the colour-magnitude and colour-colour 
diagrams, finding consistent results. 

An accurate i^'-band calibration of the red clun ip giant lum inosity was obtained and 
applied to t he red giant clump of the Galactic bulge (lAlvesI l2000l ). and of the Large Magel- 
lanic Cloud flAlves et al.ll2002t ). Our magnitudes are in the 2MASS mag nitude syste m, and 
we transformed the i^-band magnitudes of clump giants of Hipparcos (lAlved 120001 ) to Kg 
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magnitudes using -ft" = K<^ + 0-044. The resulting zero point differences should be less than 
0.02 magnitudes (lAlves et al.ll2002l ). For the Ks photometry, the distance modulus to a red 
clump giant in the outer disk would be: 



[{J-K,)-{J-K 



(1) 



where we adopted AkJ{Aj-AkJ = 0.73, {J-K,)o = 0.70±0.05 and My. = -1.65± 0.03 as 
the mean values for the red clump giants of the Milky Way disk (lAlves et al.ll2002l ). There 
should be negligible metallicity dependence of these mean values because we are looking 
at a stellar population of the Milky Way disk that should be similar to that of the Solar 
neighbourhood where Hipparcos distances of clump giants wer e calibrated. Adopti ng these 
means magnitudes and colours, and the reddening coefficients (ICardelli et al.lll989l ) yields: 



/i= -5 + 51ogrf(pc) = K, - 0.73 {J - K,) + 2.16 



(2) 




Fig. 2. — Distance distribution for clump giant stars along different lines of sight. The 
distance measured to the edge of the Milky Way stellar disk is labelled for each field. The thin 
solid line is the smoothed distribution. These stellar distributions are rich with information 
about the disk structure. For example, we discarded the presence of distant star clusters, 
where the red giant clump could mask an edge feature, but see some over-density fluctuations 
that may be due to spiral arms. 

Using this equation we computed the distance modulus (and distance in kpc) for every 
single clump giant candidate in the fields. For the UKIDSS-GPS K-hand data the com- 
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putations are the same, but with Ak/{Aj — Ak) 
Mk = -1.61 ±0.03. 



0.68, (J - K)o = 0.66 ± 0.04 and 



In order to determine the disk edge we computed the distribution in distance for the 
clump giants in each field. Thus we coadded the distance distribution of fields located above 
and below the plane at same Galactic longitude. The coadded distance distribution for 
each lon gitude was aii alyzed non-parametrically using the local likelihood density estimation 
method (lLoadeiill996l ). Finally, the derivative method was applied for the distribution curve. 
This method is simple and robust, and can be used to detect the structures in the distance 
distribution. We defined as the disk edge the outermost point of minimum derivative in each 
case. The uncertainties were calculated with a Monte Carlo procedure, assuming Poisson 
errors. 

The uncertainties due to the photometric errors do not produce any trend in the calcu- 
lations and are much smaller than the statistical error. The resulting distance distributions 
along different lines of sight are shown in Fig. |2l In all fields there is a density drop in 
the clump giant distribution at a corresponding Galactocentric distance of about 13.9 kpc 
(assuming Ro = 8 kpc), as listed in Tabled! 

We first noticed the sharp termination of the clump giant distribution in the color- 
magnitude diagram of the VVV field dOOl {l,b = 295.4°,- 



20081). 



fSaito et al. 


2010). Similar 


tic longitudes ( 


Lucas et al. 



We studied low-extinction fields across the Milky Way disk (see Tabled]). When these 
were available, we selected pairs of fields located above and below the plane at the same 
Galactic longitude, in order to account for the Galactic warp, because if the disk is warped, 
the line of sight can leave the stellar d isk before this ends . For example, at / = 300° the mean 
warp location is 1° below the plane (IRobin et al.lll992l ). so the selected VVV field d003 is 
conveniently located to probe the full extent of the disk. At most other longitudes explored 
here the warp is absent, with the mean plane being at 6 = 0° (Fig. [3]). We have also explored 
a few other UKIDSS-GPS fields located at inner longitudes that yield lower distances. These 
fields with / < 60° were discarded because they are heavily reddened and crowded, and the 
photometry is not as deep as the rest, presumably due to enhanced contamination from the 
near-side of the inner bar. Systematic errors such as variations of the reddening law have 
not been included. 



It is important to note that the present results are independent of models. They rely on 
basic assumptions such as that the Hipparcos stellar sample is representative of the entire 
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the Milky Way disk, and that clump giants are reliable distance indicators and t racers of the 



old a nd intermediate age populations. However, the Besangon Galactic model (iRobin et al. 



20031 ) with a scale-length of 2.4 kpc, including a flaring and warped outer disk, with a disk 



cutoff at i? = 14 kpc seems to reproduce well the observed distance distributions (see FigS]). 
Data and model show the same features, with a sharp increase, the peak around 5 kpc 
and the exponential decrease. The edge is detected in the data and model at the expected 
place within the errors [d = 10.1 kpc (corresponding to R = 13.9 kpc), and d = 9.5 kpc 
{R = 13.4 kpc), respectively]. 



3. Discussion and Conclusions 



A few studies have discussed the edge of the Milky Way towards the anticentre region 
previously. The old Galactic disk apparently does not extend beyond 14 kpc on the basi s 
of optical V vs. {B — V) colour-magnitude diagrams and star counts (IRobin et al.lll992l ). 
which does not contradict the presence of other more distant young stars. Models of the 
COBE-DIRBE infrared emissi on maps yield that the edge of the disk is at 4 kpc from the 
Sun (assuming Ro = 8.5 kpc, iFreudenreichI Il996l ) . DENIS star counts reveal the cut off of 
the stellar disk at a Galactocentric distance of 15 ± 2 kpc (iRuphy et al.lll996l ). 



In contrast, star counts from 2MASS reveal no radial disk truncation at 14 kpc (IMomany et al. 



20061 ) . Again, the depth of 2MASS is not enough to reach large distances, which is possible 
with the VVV and UKIDSS-GPS survey that reach 3-4 magnitudes fainter. More recently, 
it was found that the 2MASS star counts are best fit if the external disk is truncated at 



12-14 kpc (IReyle et al.l 120091 ). while early A-type stars in the anticentre from the IPHAS 
survey reveal an exponential disk out to 13 kpc, with a steeper decline beyond that distance 
(ISale et al.ll2010l ). Most of the previous evidence for cutoffs in the stellar distribution has 
been acquired at the anticentre fields. The present deep exploration of different fields in 
three galactic quadrants finds consistent results. 

Even though we interpret this termination of the clump giant distribution as a trun- 
cation, we note that a break in the slope of the distribution would also be consistent with 
the data. It is difficult to measure the sharpness of the cutoff, because the number of stars 
drop rapidly with distance, but as the surveys progress, we will have more fields in order to 
explore this issue. 

This does not mean that one could not find other young stellar sources beyond the 
distance that we measure. There are stars detected beyond the edge of the old stellar disk. 
For example, there is a population of distant young stars at i? = 20 kpc, in two fields 
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Table 1: Selected GPS and VVV fields with relatively low extinction (jSchechter et al.lll993l ) 
in separate directions across the Milky Way plane. 



Field ID 


in 


Ml 


d [kp 




R[kp 


c] 


GPS-U098+ 


98.4 


+1.0 


10.1 


± 


1.0 


13.9 


± 


1.3 


GFS-uOQO- 


90.0 


-1.7 


11.2 


± 


0.9 


13.8 


± 


1.1 


GPS-U060- 


60.0 


-2.0 


16.3 


± 


0.7 


14.1 


± 


0.6 


GPS-U060+ 


60.0 


+2.0 


16.6 


± 


0.8 


14.3 


± 


0.7 


VVV-d028 


334.7 


-1.9 


20.3 


± 


0.7 


13.5 


± 


0.5 


VVV-dl42 


334.8 


+ 1.9 


20.7 


± 


1.0 


13.8 


± 


0.7 


VVV-d014 


314.3 


-1.9 


17.4 


± 


0.9 


13.1 


± 


0.6 


VVV-dl28 


314.3 


+1.8 


16.7 


± 


1.3 


12.6 


± 


1.0 


VVV-d003 


298.3 


-1.9 


16.4 


± 


0.8 


14.6 


± 


0.8 


VVV-dll5 


295.3 


+1.9 


16.3 


± 


0.8 


14.5 


± 


0.7 



located in the thir d quadrant, slightly above the plane (at 6 = 7° and b = 4°), based on 
UBV photometry (ICarraro et al.ll2010l ). We see no evidence of a large population of clump 
giants at that distance in our fields. 



i^'-band light traces mass in the disks of spiral galaxies (iRix fc Riekd Il993l ). and the 
integrated i^'-band luminosity is dominated by red giants. Thus the clump giants are ideal 
tracers to define overall features and structural parameters of our Milky Way, like the edge 
of its disk. The red clump giants trace the old and intermediate age populations, which 
are more uniformly distributed than young stars, following the mass of the galactic disk. 
However, thick disk and halo giants are expected to be loc ated beyond the edg e of the stellar 
thin disk, as have been detected by the SEGUE project ( Ide Jong et al.l 120101 ) . contributing 
to the background that we see in all fields beyond a Galactocentric distance of 13.9 kpc. A 
similar result has been recently obtained independently by the GLIMPSE team (Benjamin 
et al. 2011, in preparation). 

Finally, we are able to determine not only the extension but also the shape of the stellar 
disk of our galaxy for the first tirne. As ymmetric features such as lopsidedness are common 
in spiral galaxies (IRix fc Zaritskylll995l ). However, the results shown in Fig. E] suggest that 
the stellar disk of the Milky Way is not significantly lopsided. 
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Fig. 3. — Measured positions of the edge of the Milky Way stellar disk. Top panel: 
Schematic face-ori view of the Milky Way disk based on radio, infrared, and visible data 
(IChurchwell et al.ll2009( ). showing the positions of the edge of the stellar disk measured in 
different directions. The fields to the left are from the UKIDSS-GPS, the fields to the right 
are from the WV survey, and the anticentre point at i? = 14 kpc is from (IRobin et al.ll2003l ; 
Reyle et al.ll2009l ). We adopted a distance to the galactic centre of -Ro = 8 kpc. This figure 



shows that a disk radius of 13.9 ± 0.6 kpc (dotted circle) fits the data very well. Bottom 
panel: Edge-on view of the Milky Way in Galac tic /, h coordinates w ith the position of the 
warped stellar disk marked with white squares ( iMomany et al.ll2006l ). The positions of our 
fields are marked with triangles and circles, and the vertical scale has been stretched for 
clarity. 
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Fig. 4. — Comparison between data and model. Black solid line: Distance distribution 
for the field u098+ using UKIDSS-GPS data. Dashed red line: Simulation of the same 
field u098+ using the Besangon Galactic model with a scale length of 2.6 kpc and an edge 
at 14 kpc. The model fits all features seen in the data distribution and the disk edge is 
detected at the expected place within the errors (see text). The distributions are arbitrarily 
normalized, and at shorter distances the discrepancy is due to photometric saturation. 



